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ABSTRACT 

The deuterium abundances inferred from observations of the interstellar medium within 1-2 
kpc of the Sun range over a factor of three and the corresponding oxygen abundances show 
an even larger dispersion. While the lower D (and O) abundances likely result from deple- 
tion onto dust, the higher D abundances are consistent with the BBN-predicted primordial D 
abundance and chemical evolution models of the Galaxy with infall of primordial or nearly 
primordial material. The large ranges in deuterium and oxygen abundances suggest that the 
effects of depletion and/or infall have not been homogenized in the local interstellar medium. 

Key words: Galaxy: evolution - ISM: abundances. 



1 INTRODUCTION 

Of the light nuclides produced in astrophysically interesting abun- 
dances (D, 3 He, 4 He, 7 Li) during Big Bang Nucleosynthesis 
(BBN), deuterium is the baryometer of choice for several reasons. 
For example, the BBN-predicted primordial deuterium abundance 
is sensitive to the baryon (nucleon) density [yop = 10 5 (D/H)p oc 
rj^g' 6 where r/io = 10 10 (riB/n 7 )o] and the post-BBN evolution 
of D is simple and monotonic: as gas cycles through stars, deu- 
terium is only destroyed (Epstein, Lattimer & Schramm 1976; Pro- 
danovic & Fields 2003). A consequence of this simple evolution is 
that the post-BBN deuterium abundance is a probe of the fraction 
of the interstellar gas which has never been cycled through stars 
(Steigman & Tosi 1992, 1995). Models of the chemical evolution 
of the Galaxy, constrained by observations of the stellar and gas 
phase metal abundances, follow the history of the gas. As a bonus 
they predict the time-evolution of the deuterium abundance. Until 
relatively recently observations of deuterium in the Universe were 
limited to the interstellar medium (ISM) of the Galaxy, so that de- 
terminations of primordial deuterium required the input of chemi- 
cal evolution models and, hence, were model-dependent (Steigman 
& Tosi 1992, 1995; Vangioni-Flam, Olive & Prantzos 1994; Prant- 
zos 1996; Tosi 1996; Scully et al. 1997; Tosi et al. 1998). This 
changed with the advent of large, ground-based telescopes and with 
the Hubble Space Telescope (HST), which have enabled determina- 
tions of the D abundance along the lines-of-sight (LOS) to a few (at 
present, six) high-redshift, low-metallicity QSO Absorption Line 
Systems (QSOALS; Kirkman et al. 2003; O'Meara et al. 2006). At 
the same time observations of the angular distribution of the tem- 
perature fluctuations in the Cosmic Background Radiation (CBR), 
in combination with new data on the large-scale structure of the 
Universe (see, e.g., Spergel et al. 2006; Tegmark et al. 2006), have 
constrained the baryon abundance to ~ 3 %, rjio = 6.08 ± 0.19, 
leading to a BBN-predicted primordial D abundance accurate to 
~ 5 - 6 %, yop = 2.61 ± 0.15. Starting with this primordial abun- 



dance, what do Galactic chemical evolution (GCE) models predict 
for the current local ISM deuterium abundance and how well do 
these predictions compare to observations? 

This paper is organized as follows. In Section 2 the predicted 
and observed abundances of D and O in the local ISM are reviewed 
and compared, accounting for depletion onto dust grains as well as 
localized infall of metal-poor, D-enhanced gas. The implications of 
ISM inhomogeneities are further discussed in Section 3. Finally, in 
Section 4 we draw our conclusions. 



2 COMPARING ISM PREDICTED AND OBSERVED 
ABUNDANCES 

As gas cycles through stars, D is destroyed and the abundances 
of the "metals", in particular oxygen, increase (from zero primor- 
dially). In addition to star formation, gas flows play a fundamental 
role in determining the chemical composition of gas and stars in 
galaxies. In particular, infall of primordial, or nearly primordial, gas 
plays a crucial role in the chemical evolution of our own Galaxy, 
enhancing the ISM D abundance, restoring it closer to its primor- 
dial value (see, e.g., Tosi 1988a, b). As a result the ISM deuterium 
abundance is not expected to depart from its primordial value until 
the ISM metallicity increases to <; 0.1 solar (Romano et al. 2006). 
Underlying GCE models is the "conventional wisdom" that the lo- 
cal ISM is chemically homogeneous (i.e. the timescale for mixing 
in the local ISM is much shorter than stellar/chemical evolution 
timescales). Consequently, GCE models predict mean abundances 
for deuterium and the metals, e.g. oxygen, in the vicinity of the 
Sun at the present time. The data obtained from the Far Ultraviolet 
Satellite Explorer (FUSE; Oliveira et al. 2006; Oliveira & Hebrard 
2006), however, reveal a factor three range for the D abundance and 
an even larger range for the O abundance, clearly in conflict with 
this simplifying assumption. The gas phase abundances of D and O 
are currently not distributed homogeneously in the local ISM. 
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2.1 Predicted ISM abundances 

The GCE model-predicted ISM abundance of D is related to the 
primordial D abundance through the "astration" factor, /d = 
2/dp/2/d M - Many GCE models are found in the literature which 
address D evolution in the Milky Way. Generally, they adopt dif- 
ferent prescriptions for the star formation rate, infall and outflow 
rates, initial mass function (IMF) and stellar lifetimes (see, e.g., 
Vangioni-Flam et al. 1994; Prantzos 1996; Tosi 1996; Scully et 
al. 1997; Tosi et al. 1998). Successful models which account for 
a wealth of observational data always find "low" values for /b , in 
the range ~ 1.5 — 3. Recent models (Romano et al. 2006) in which 
the prescriptions of the stellar IMF and stellar lifetimes (Scalo 
1986; Kroupa, Tout & Gilmore 1993; Maeder & Meynet 1989; 
Schaller et al. 1992) are allowed to change, predict that the astration 
factor and the current ISM oxygen abundances lie in the ranges, 
1.39 <> fo <, 1.83 and 7 <> yo £ 12, where yo = 10 4 (O/H). 
For a BBN-predicted primordial D abundance of yop = 2.61, the 
GCE-predicted ISM D abundance is expected to lie in the range, 
1.43 <S yE M £ 1.88. 

There are at least two physical effects capable of accounting 
for the observed variations in D and O abundances along different 
LOS within ~ 1 — 2 kpc of the Sun. For an explanation of the vari- 
able deuterium abundances Linsky et al. (2006) embrace the sug- 
gestion of Draine (2004), based on earlier work of Jura (1982), that 
deuterium may be preferentially depleted (compared to hydrogen) 
onto dust grains. While Linsky et al. do not address the variation in 
oxygen abundances, depletion might provide a viable explanation 
for the dispersion among the observed oxygen abundances as well. 

Another possibility is recent infall to the disk of the Galaxy 
of unprocessed (D-rich, O-poor) material which has been incom- 
pletely mixed in the ISM (see also Geiss, Gloeckler & Charbonnel 
2002). In this case there would be an anti-correlation between the 
D and O abundances observed along such contaminated LOS. For 
example, if yjf M and y^ M = i/ DP //d are the "true" ISM abun- 
dances and 2/§ BS an d 2/d BS the abundances observed along a spe- 
cific LOS, then incompletely mixed infall predicts 



OBS 
2/D 



. / ISM \ una i is 

2/dp + (2/d -yDPjyo IVo 



OBS / ISM 



(1) 



For example, if the Scalo (1986) IMF along with the stellar life- 
times of Schaller et al. (1992) are adopted, then according to Ro- 
mano et al. (2006), / D = 1.39 and y§ M = 8.4, so that 
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is predicted. 

Since FUSE is often unable to provide sufficiently accurate 
H i column densities, to interpret these data people usually rely on 
H i column density determinations from other telescopes/observers. 
As a result, it is possible that some of the dispersion among the 
observed D and O abundances may be due to systematic errors in 
the associated Hi column densities (Timmes et al. 1997; Hebrard 
& Moos 2003; Steigman 2003). However, the ratio of D to O is 
unaffected by such errors. For 10(D/O) = yn/yo, eq. 2 predicts 

10(D/O) = 26.1/j/g BS - 0.87. (3) 

Of course, it may well be the case that both mechanisms (de- 
pletion and infall) are at work in the ISM. If so, the anticorrela- 
tions described by equations 1 — 3 provide an upper envelope to 
the observations in the j/d — yo and D/O - yo planes (see Figs. 3 
and 4). Infall of (nearly) primordial gas is often invoked to account 
for abundance variations in the ISM (e.g. Cartledge et al. 2004; 
Knauth, Meyer & Lauroesch 2006, and references therein). The 



required infall episodes should have occurred recently enough to 
avoid complete mixing of the ISM. The exact timescale on which 
clouds of partially ionized gas falling into the magnetized ISM will 
be thoroughly homogenized is not clear. At the same time it is clear 
that if depletion is the correct explanation for the observed varia- 
tions in D and O abundances, then local depletions have not been 
homogenized either. Therefore, although it will be seen below that 
recent inhomogeneous infall on small (pc-sized) scales is not re- 
quired to achieve concordance between the observed ISM abun- 
dances and GCE models, neither should it be rejected out of hand 
as contributing to the observed dispersion among the gas phase D 
and O abundances. 



2.2 Observed - gas phase - ISM abundances 

The abundance of deuterium in the relatively local ISM (within 1 — 
2 kpc of the Sun) has been explored for some 30 years beginning 
with observations by the Copernicus satellite and continuing with 
a series of space-based UV detectors (e.g. IUE, GHRS and STIS on 
the HST) up to the current IMAPS instrument on the FUSE satellite. 
The most recent and comprehensive abundance analyses are based 
on the data acquired by FUSE (Oliveira et al. 2006; Linsky et al. 
2006; Oliveira & Hebrard 2006). These data are adopted in our 
analysis here, together with the abundances in the warm neutral 
medium of the lower Galactic halo recently determined by Savage 
et al. (2007). 

Linsky et al. (2006) suggest that local variations in the ob- 
served gas-phase D/H ratios can be explained by spatial variations 
in the depletion of deuterium onto dust grains. While this hypothe- 
sis of localized depletion provides a reasonable explanation for the 
observed variations in the D and O abundances, there is little direct 
evidence in support of it. For example, although the FUSE obser- 
vations cover a range of four orders of magnitude in Hi column 
density [17.6 <5 logN(H i) <J 21.2], both the lowest and the high- 
est deuterium abundances are found for 19.2 <^ logN(H i) <, 21.0 
(see Oliveira et al. 2006, their fig. 17, upper panel, and Oliveira 
& Hebrard 2006). Of course, since the LOS explored by FUSE 
range in distance (to the background star) from tens of parsecs to 
1 — 2 kpc, high H i column density is not equivalent to high gas 
density. With few exceptions the FUSE data do not permit a direct 
determination of the local gas density. A surrogate for the gas den- 
sity, n(H i)=N(H i)ld, where d is the distance to the background il- 
luminating star, ranges over nearly three orders of magnitude, from 
~ 0.002 cm~ 3 to ~ 0.6 cm~ 3 . While this estimate doesn't pre- 
clude much higher local densities, the largest of these is smaller 
than the densities associated with molecular clouds where deple- 
tion is expected to be efficient. Furthermore, here too, the evidence 
in support of D-depletion is ambiguous in that although some of 
the lowest (depleted?) D abundances are found at the highest val- 
ues of n(H i), so too are the majority of the highest (undepleted?) 
D abundances (see Oliveira et al. 2006, their fig. 19, upper panel). 
Alternatively, if dust depletion is, indeed, the culprit it might be 
anticipated that the deuterium and oxygen abundances should be 
anticorrelated with the LOS reddening, as measured by E(B— V) 
or, with the molecular hydrogen fraction. We explored these possi- 
bilities (see Fig. 1 for deuterium) using estimates of E(B— V) from 
Schlegel, Finkbeiner & Davis (1998) and H2 data from Hoopes et 
al. (2003), Linsky et al. (2006), Oliveira & Hebrard (2006). We 
find that although there is a tendency for the lowest (highest) deu- 
terium abundances to correspond to the highest (lowest) E(B— V) 
values and H2 fractions, many of the least reddened LOS, as well as 
those with the smallest H2 fractions, where depletion may be less 
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E(B-V) log N(H a ) 

Figure 1. D/H versus E(B— V) (left panel) and logN(H2) (right panel). Esti- 
mates of E(B— V) are from Schlegel et al. (1998); H2 data are from Hoopes 
et al. (2003), Linsky et al. (2006) and Oliveira & Hebrard (2006); D/H deter- 
minations are taken from Oliveira et al. (2006; dots) and Oliveira & Hebrard 
(2006; squares). 



likely, also have low D abundances. Furthermore, two out of eight 
LOS with the highest H 2 fractions (HD41161 and PG 0038+199) 
have some of the highest D abundances observed. Similar trends 
are found for the oxygen abundances [low O/H at low E(B— V), 
low H2 fraction]. 

Perhaps the strongest evidence in support of D deple- 
tion comes from the correlation between D and Ti abundances 
(Prochaska, Tripp & Houk 2005; but see also Hebrard et al. 2005). 
The FUSE data (Oliveira et al. 2006; Linsky et al. 2006; Oliveira 
& Hebrard 2006) reveal a similar trend between D and Fe. These 
data, plotted in Fig. 2, show that the lowest D abundances are, in- 
deed, correlated with the lowest Fe abundances [log yp c ^ 0.0, 
where i/f c = 10 6 (Fe/H)]. However, these data also suggest that 
deuterium may be undepleted along several of the LOS for which 
log yp c 0.0 — 0.2. This calls into question the Linsky et al. 
(2006) assumption that a lower limit to the "true" ISM D abun- 
dance is provided by the maximum of y§ BS . 

If dust depletion is the dominant mechanism responsible for 
the observed dispersion among the gas-phase abundances of deu- 
terium and oxygen, the D and O abundances should be correlated 
(as D appears to be with Ti and Fe). Incompletely mixed infall 
of unprocessed material would give rise to an anticorrelation be- 
tween the D and O abundances. To explore these trends in the 
FUSE data we concentrate on those LOS with data (mostly from 
Oliveira et al. 2006) for both the D and O abundances. The D/H 
and D/O ratios versus O/H are plotted in Figs. 3 and 4. If deple- 
tion dominates, none of the observed abundances should exceed the 
GCE-predicted values. For the fiducial model we adopted the Scalo 
(1986) IMF along with the Schaller et al. (1992) stellar lifetimes, 
so that / D = 1.39 and t/o SM = 8.4, leading to yff M = j/dp//d = 
1.88 ± 0.11. The upper panel of Fig. 3 shows D/H versus O/H 
along with the upper bounds to their GCE-predicted ISM abun- 
dances if depletion dominates. We note that while most of the LOS 
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Figure 2. ISM deuterium abundances versus iron abundances yp c = 
10 6 (Fe/H). Data from: Linsky et al. (2006) for Fe and Oliveira et al. (2006) 
for D (dots); Oliveira & Hebrard (2006; full squares); Savage et al. (2007; 
upside-down triangle). The star shows the abundances predicted by the 
fiducial chemical evolution model (see the text). D and Fe depletion re- 
quire that the data lie below and to the left of the dashed lines. The solid 
curve shows the effect of local abundance fluctuations resulting from in- 
completely mixed infall of primordial (D-enhanced, Fe-free) material. 

have D or O abundances which lie below these bounds, consistent 
with depletion, none of them are more than la above these GCE- 
predicted ISM abundances, consistent with some of these LOS be- 
ing depletion-free. One possible exception may be the LOS towards 
LSE44 (Friedman et al. 2006; see Fig. 3). However, as Liszt (2006) 
notes, it is important to account for molecular hydrogen when de- 
termining the ISM deuterium abundance, since failure to do so bi- 
ases the resulting D/H (and O/H) ratio(s) upward. Liszt notes that 
molecular hydrogen is unlikely to be negligible (compared to H 1) 
whenever N(H 1) > 4 x 10 20 cm -2 and should be considered when- 
ever N(Hi) > 1.4 x 10 20 cm -2 . This makes the LOS to LSE44, 
with N(H 1)= 3.3 x 10 20 cm~ 2 , a good candidate for possibly over- 
estimated D and O abundances*. This possibility is indicated by 
the arrows in Figs. 3 and 4; note that the D/O ratios are unaffected. 

The large difference in H 1 and D 1 column densities, ~ 5 or- 
ders of magnitude, is a potential source of systematic errors for 
the D/H determinations. The difficulties associated with determin- 
ing D/H may be avoided by measuring the D/O ratio. Since charge 
exchange reactions between singly-ionized and neutral H and O 
(Field & Steigman 1971) couple Hi to Oi, oxygen (Oi) is a good 
tracer of hydrogen (H 1), so that the D/O ratio might be a good 
proxy for the D/H ratio (Timmes et al. 1997; Hebrard & Moos 
2003; Steigman 2003). The D/O ratio is plotted (versus O/H) in 

* Although within the errors the D and O abundances measured towards 
LSE44 are, in fact, consistent with the predictions of our fiducial GCE 
model (j/p M = 1.88, j/q M = 8.4), it is worth emphasizing that if molec- 
ular hydrogen were to contribute N(H2)/N(H 1) = 0.19 along this LOS, then 
j/D would be reduced from 2.24 to 1.88, and yo from 11.3 to 9.5, in even 
better agreement with the model predictions. 
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Figure 3. The upper panel shows the deuterium versus oxygen abundances. 
The lower panel plots the ratio of deuterium to oxygen as a function of 
the oxygen abundances. In both panels, the FUSE data are displayed as 
filled circles (Oliveira et al. 2006), squares (Oliveira & Hebrard 2006) and 
the upside-down triangle (Savage et al. 2007), while the star indicates the 
ISM abundances predicted by the fiducial chemical evolution model (see 
the text). The horizontal dashed and dotted lines in the upper panel are the 
predicted ISM deuterium abundance and its ±1<t uncertainty, respectively. 
The vertical dashed lines in both panels correspond to the predicted ISM 
oxygen abundance. The dashed curve in the lower panel shows the maxi- 
mum values of D/O if O is depleted, but D is not. Data below and to the 
left of the dashed lines are consistent with D and O depletion. The arrows 
in each panel indicate the effect of a possible contribution of molecular hy- 
drogen to the D and O abundances for the LOS towards LSE 44. 

the lower panel of Fig. 3. Indeed, the FUSE data (Oliveira et al. 
2006; Linsky et al. 2006; Oliveira & Hebrard 2006) reveal that the 
D/O ratio has an even larger dispersion than that for D/H (a factor 
of ~ 5, similar to that for O/H), suggesting that the variations in O 
and in D abundances may be real and not the result of systematic 
errors. However, note that, as seen in the lower panel of Fig. 3, the 
upper envelope of the D/O versus O/H data is consistent with an 
undepleted ISM D abundance of y§ M = 1.88. 

While the lower D and O abundances shown in Figs. 3 and 4 
may have resulted from the depletion of gas phase D and O, the 
highest deuterium and oxygen abundances derived from the FUSE 
data are consistent with present epoch average ISM deuterium and 
oxygen abundances of y^ M = 1.88 ± 0.11 and y^ M = 8.4. For 
deuterium, these are consistent with a primordial abundance yup = 
2.61 ± 0.15 and a Galactic astration factor of 1.39, suggesting that 
gas phase D (and O) may be undepleted along at least some LOS 
in the local ISM. Modest depletion of deuterium along some of 
these high-D LOS could also be consistent with the data provided 
that the D (and O) variations are due, in part, to the contribution 
of incompletely mixed, D-enhanced, O-reduced infall^ . This latter 
possibility is shown in the two panels of Fig. 4. Even including the 

T Notice that our models are computed assuming infall of primordial (i.e. 
O-free) material. The conclusions are unchanged if the infalling material is 



Figure 4. Deuterium abundances and D/O ratios as in Fig. 3. The solid 
curves show the effect of local fluctuations in the abundances in the absence 
of depletion but accounting for infall of primordial (D-enhanced, O-free) 
gas on small spatial scales. 

LOS to LSE 44, whose D/H and O/H ratios might be overestimated 
if there is a non-negligible contribution from unobserved H2, the 
FUSE data are consistent with primordial deuterium, GCE models, 
depletion, and - possibly, but not necessarily - local fluctuations 
due to infall. 



3 DISCUSSION 

Despite the evidence of 30 years of data to the contrary, until re- 
cently the conventional wisdom has been that the local ISM is well 
mixed, with unique abundances of the chemical elements, D and O 
in particular. The FUSE data, complementing earlier observations, 
which reveals a factor of three range for D/H (0.7 yo 2.2) 
and an even larger dispersion among the O/H ratios (2 ^ yo ^ 
8—11), lays to rest this overly simplified assumption. While Linsky 
et al. (2006) don't address the spread in oxygen abundances, they 
invoke localized depletion to account for the observed range of deu- 
terium abundances. Although independent observational evidence 
for depletion is currently ambiguous, this mechanism does provide 
a description of the data. If so, the question becomes which, if any, 
of the observed D/H ratios represents the "true" (undepleted) ISM 
abundance? 

Linsky et al. (2006) assume that the true D abundance is 
"equal to or slightly above the highest measured (D/H) gas ratios". 
Limiting themselves to the 5 LOS with the highest D abundances, 
Linsky et al. conclude that y§ M ^ 2.17 ± 0.17 (or, including a 
correction for the Local Bubble D abundance, y§ M > 2.37 ± 0.24 
or yjj = 2.31 ± 0.24, depending on the scheme for weighting the 
data). These choices correspond to an astration factor ~ 1.1 — 1.2, 
considerably smaller than those predicted by GCE models (Tosi 

slightly enriched (^i n f a ii $5 0.2 Zq, as for instance that coming from the 
Magellanic Stream; Tosi 1988b; Tosi et al. 1998). 
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1996 and references therein; Chiappini, Renda & Matteucci 2002; 
Romano et al. 2006). However, given the errors associated with the 
FUSE D/H determinations, it is unclear why their estimate should 
have been limited to the 5 LOS with the highest central values of 
the D/H ratios since there are more than a dozen other LOS whose 
central values of D/H are within 1 a of them. To test our hypoth- 
esis that an astration factor of 1.39, in combination with a BBN 
expected primordial abundance of j/dp = 2.61 ± 0.11, is consis- 
tent with the FUSE data, we have considered the 18 FUSE LOS 
with the highest D/H ratios (j/d ^ 1.5). Of these, nine have central 
values yo ^ 1.9 and nine have yo ^ 1.7. Using these data we may 
test the assumption that the undepleted ISM D abundance is con- 
sistent with y rj SM = 1.88. We find that 13/18 (72%) have central 
values of D/H within 1 a of our adopted ISM abundance and none 
are more than 1 .4 <r away. The reduced X 2 is 0.70 (for 17 de grees 
of freedom) suggesting a more than ~ 70% probability that the D 
abundances along these 18 LOS are consistent with our adopted 
value (yo M = 1.88). 

Indeed, we have been "conservative" here in that, had we 
adopted for the primordial D abundance the O'Meara et al. (2006) 
estimate of j/dp = 2.84 inferred from the data for 6 high-z, low-Z 
QSOALS, consistency with the FUSE data is improved, allowing a 
slightly larger astration factor of fo ~ 1.5. 

In addition to inhomogeneously distributed depletion and pos- 
sible systematic errors in some of the H i column density determi- 
nations, some of the variations in D and O abundances could be 
due to incompletely mixed, recent infall of nearly primordial gas. 
As shown in the lower panel of Fig. 4, the upper envelope of the 
D/O ratios reveals a trend with O/H consistent with localized in- 
fall. However, although infall may contribute to the large disper- 
sions in the D and O abundances, depletion appears to be required 
to account for the lowest D/H and O/H observed ratios. In any case, 
with or without localized infall, at least some of the current GCE 
models are consistent with the FUSE-inferred and CBR-predicted 
BBN D abundances. 



4 SUMMARY AND CONCLUSIONS 

Standard galactic chemical evolution models predict the climate 
and not the weather and, hence, they provide average abundances 
for deuterium and oxygen (as well as other chemical species) in 
the local ISM. This is in contrast to the large dispersion among 
the gas phase D and O abundances in the neutral ISM revealed 
by 30 years of data. Here we have demonstrated, despite recent 
claims to the contrary (Linsky et al. 2006; Savage et al. 2007), that 
current models for the chemical evolution of the Milky Way (Ro- 
mano et al. 2006) are consistent with the FUSE inferred and CBR- 
predicted BBN deuterium abundances. While depletion of D (and 
O) onto dust grains (Draine 2004) alone may suffice to account 
for the observed variations, recent infall of (nearly) primordial gas 
on small spatial scales can not be ruled out as a viable explana- 
tion for at least some of the highest observed D abundances. It is 
worth emphasizing that such inhomogeneous infall - which must 
have occurred recently enough that mixing processes have not yet 
erased its elemental abundance signature - sums up to the continu- 
ous, homogeneous infall required by chemical evolution models to 
account for the solar neighbourhood metallicity data. Finally, no- 
tice that neglect of possible contributions from molecular hydrogen 
(Liszt 2006) may have biased some of the observed abundances up- 
wards, enlarging the observationally inferred range for the present- 
day D/H (and O/H). 
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